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Abstract 

We obtain a consistency relation for the observed three-point correlator of galaxies. It in¬ 
cludes relativistic effects and it is valid in the squeezed limit. Furthermore, the consistency 
relation is non-perturbative and can be used at arbitrarily small scales for the short modes. 
Our results are also useful to compute the non-linear relativistic corrections which induce a 
signal in the observations that might be misinterpreted as primordial non-Gaussianity with a 
local shape. We estimate the effective local non-Gaussianity parameter from the relativistic 
corrections. The exact value depends on the redshift and the magnification bias. At redshift 
of z = 1.5, in the absence of magnification bias, we get ~ —3.0. 


1 Introduction 

The experimental developments of large-scale galaxy surveys over the past few decades, in 
addition to high precision maps of the Gosmic Microwave Background (CMB) have improved 
extensively our understanding of the universe over a large range of scales. With the upcom¬ 
ing galaxy redshift surveys, such as Euclid [1], the Dark Energy Spectroscopic Instrument 
(DESI) [2], the Large Synoptic Survey Telescope (LSST) [3] and the Wide-Eield InfraRed 
Survey Telescope (WEIRST) [4], covering progressively a larger fraction of the sky and deeper 
redshift range, this trend will continue further in the future. 

The enhanced precision of these surveys requires high accuracy theoretical predictions 
for the observed quantities rather than theoretically convenient quantities. Moreover, on the 
very large-scales probed by these surveys, the Newtonian description of galaxy clustering is 
not valid. The relativistic description of the observed galaxy number density, in terms of the 
observed redshift and angle of galaxies, have been developed up to first order in perturbation 
theory [5-9], which has allowed them to compute the impact of relativistic effects on the 
galaxy power spectrum. The galaxy bispectrum contains much additional information not 
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present in the power spectrum as there are many more modes and it suffers from different 
systematics. The computation of this bispectrum in principle needs the use of second order 
perturbation theory which is very involved and non-trivial. Unlike the CMB, where the 
non-linearities are small and therefore (for Gaussian initial conditions) the fluctuations are 
nearly Gaussian, the large-scale structure statistics is intrinsically non-Gaussian due to the 
gravitational collapse. One needs therefore to go beyond the two-point statistics and hence 
beyond the linear level in perturbation theory. Recently, several groups have performed the 
second-order calculation [10-13]. 

The expression for the observed galaxy three-point correlator (bispectrum) obtained 
from the second-order calculation is rather complicated and has contributions from various 
terms. Moreover the region of validity of this calculation is only in the weakly non-linear 
regime. 

In this paper we point out that a consistency relation for the observed galaxy bispectrum 
can be obtained in a non-perturbative manner in the squeezed limit (that is one wavenumber 
much smaller than the other two in Fourier space). This consistency relation allows to express 
the observed bispectrum in terms of the observed galaxy power spectrum. This is done by 
correlating the observed galaxy power spectrum in the presence of a long mode with the 
observed galaxy overdensity on that large wavelength. This is exactly computable since the 
effect of a long mode on the short-scale dynamics is simply a change of frame (or a residual 
gauge transformation). 

The consistency relations are a consequence of the fact that the system at hand pos¬ 
sesses non-linearly realized symmetries. A familiar example of linearly realized symmetry 
is the translation invariance of n-point correlation functions of, for example, the dark mat¬ 
ter overdensity field 5. Non-linearly realized symmetries, such as time dependent spatial 
transformations, lead to more interesting consequences. The first example of a non-trivial 
consistency relation resulting from a non-linearly realized symmetry was pointed out by Mal- 
dacena in the context of single field inflation [14] (see also [15]). In this case the consistency 
relation is a consequence of the spatial dilation which is non-linearly realized by the curva¬ 
ture perturbations C- More recently, inflationary consistency relations for other non-linearly 
realized symmetries have been derived [16-21]. 

In the context of the large-scale structure, it was first pointed out by Kehagias & Riotto 
[22] and Peloso &: Pietroni [23] that there are non-linearly realized symmetries enjoyed by 
the system of equations (the continuity, the Euler and the Poisson equations), governing a 
pressureless fluid coupled to gravity in the Newtonian limit. The key point is that there 
are large-scale structure observables, e.g. the dark matter peculiar velocity, which transform 
non-linearly under some symmetry transformation, while the same transformation shifts the 
dark matter overdensity linearly. This generates consistency relations where the peculiar 
velocity, and not the overdensity, plays an analogous role to that of the soft-pion in the 
well-known soft-pion theorem. The consistency relations are in general between the squeezed 
limit {n + l)-point function and the n-point function of dark matter/galaxy overdensities. In 
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the general relativistic case, these non-linearly realized symmetries that involve the change of 
the coordinates can be understood as the result of a residual gauge symmetry. The extension 
of the consistency relations to the relativistic case has been done in Refs. [24, 25]. 

The advantage of using the consistency relations in the large-scale structure is that they 
are directly valid for the observed (and therefore gauge-invariant) quantities, in the same way 
they are for the CMB bispectrum [26-29], and this is what we are going to exploit in the rest 
of the paper. Our results are useful as a consistency check for the analytical second-order 
computations of the bispectrum of the observed galaxy overdensity and provide a simple and 
easy way to produce the squeezed limit. 

Furthermore, since some of these relativistic effects in galaxy number counts give rise 
to the three point function that mimics that of a local primordial non-Gaussianity (for a 
review, see Ref. [30]), removing their contribution is important for obtaining constraint on 
primordial non-Gaussianity from galaxy surveys. We will estimate the effective local non- 
Gaussianity parameter from the non-linear General Relativity (GR) corrections. The exact 
value depends on the redshift and the magnification bias, see figure 2. At redshift of z = 1.5 
and in the absence of the magnification bias we get ~ —3.0. 

This paper is organized as follows. In section 2 we review the coordinate transformation 
that induces the long wavelength mode and continue in section 3 to describe the effect of the 
long wavelength perturbation mode on the observed galaxy overdensity. We use this result 
to calculate the bispectrum in the squeezed limit by correlating two short-mode overdensities 
with a long wavelength mode. We then decompose the bispectrum in terms of spherical 
harmonics in section 4. In section 5 we estimate the value of the effective induced by 
these relativistic corrections. Finally, we conclude in section 6. 

2 Adiabatic modes and residual gauge symmetry 

The squeezed limit of an (n-|- l)-point correlation function in which one of the scales is much 
larger than the other n can be written as the correlation of a long-wavelength mode with the 
n-point function in the presence of that long mode, e.g. 

(cI)(Ai)Ag(Ai)... Ag(AO)^^^^^^ = ($(Ai)(Ag(Ai)... Ag(A„)) J , (2.1) 

where 4> is Bardeen’s potential, Ag is the galaxy number overdensity, A is the characteristic 
scale of each variable, and (• • • )l is to be interpreted in the conditional probability sense: 
it is an average given the condition that when smoothed at large-scales Bardeen’s potential 
takes the value ‘^(Al). 

It thus suffices to compute the effect of the long-wavelength potential on the short scale 
dynamics. This task is greatly simplified by the fact that on large enough scales, one can 
approximate Bardeen’s potential as a Taylor expansion 

4>L(r, x) = 4>L(r, xs) + {x - xs) ■ V4>L(r, xs) , (2.2) 
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where we have centered the Taylor expansion around an arbitrary point close to where we 
evaluate the short-scale n-point function. We will now show that in a ACDM universe this 
large-scale potential is a gauge mode: an appropriate coordinate transformation can set these 
first two terms of the Taylor expansion to zero. We will summarize the derivation given in 
Refs. [16, 24]. 

We will work to first order in the gradient of the long mode, which means we will ignore 
its curvature and tidal forces didj^L- Though these effects are in general large at 

sub-horizon scales, their effect would not be confused in observations with a violation of the 
consistency relation connecting the squeezed limit of the 3-point function with the power 
spectrum. For simplicity, we also choose to be perturbative in the short-mode potential 
and velocities ~ di^s which are expected to be small even at very small scales. Note 
that we never assume the short-mode density ~ V^4>£- to be small and in this sense we are 
non-perturbative. 

Let us start by writing the perturbed FLRW metric in Poisson gauge 


ds^ = a^(r) 


(1 -h 24>)dr^ + (1 - 2T)df^ , 


(2.3) 


where and T are Bardeen’s potentials and the Poisson gauge condition fixes the time-space 
part of the metric to be zero and the tensor part (which we ignore for simplicity) to satisfy 
didj'y'^^ = 0. Now, consider the following coordinate transformation 

T T = T + e(r) -|- X • (r) -|- a{T, x) , (2-4) 

X* X* = x*(l -I- A -|- 2x • 6) — 6*x^ -|- C*(r) (2.5) 

where a satisfies = —2di{^ + 4')^'*, A and b are constants, and e and ^ are functions of 
time. 

The key point is that, after performing this transformation, the metric is still of the 
form given by Eq. (2.3). The potentials are now 

^ = + x-i"+ n{e + x-^)], (2.6) 

^ = -[X + 2x-b + n{e + x-^)]. (2.7) 


Note that the transformation leaves the dynamics unchanged, while it changes perturbation 
quantities. Therefore, since these perturbations are still described by the same gauge condi¬ 
tion, the symmetry at issue is a non-linearly realized symmetry. The transformation we have 
performed is a residual gauge symmetry of the Poisson gauge^. This means that even when 
one completely fixes the gauge to be the Poisson gauge in which the quantities appearing in 
the metric are Bardeen’s potentials, there is a residual gauge freedom under some coordinate 
transformations that do not leave the background invariant, such as Eqs. (2.4, 2.5). Erom 
now on will ignore o. in Eq. (2.4) since its effect is of higher order under our approximations. 

residual gauge symmetry is of course there for any other gauge different from the Poisson one, which we 
have chosen because it is usually adopted as a starting point in the computation for the large-scale structure. 
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An arbitrary coordinate transformation that does not leave the background FLRW 
metric invariant can be interpreted as inducing perturbations of the metric which would have 
no physical meaning whatsoever. However, in a context where the curvature perturbation ( 
is conserved, such as for ACDM even deep inside the horizon, the transformations (2.4, 2.5) 
can be chosen to solve the large wavelength limit of Einstein’s equations. In order to satisfy 
this, the transformation parameters obey the following conditions [16, 24], 

VL = -(e + x-^), ( 2 . 8 ) 

where v is the velocity potential u* = diV. We see that corresponds to a long-wavelength 
velocity mode. Additionally we require, 

e' + 2ne = -X , + 2ni = -2b. (2.9) 

In order to see when these conditions are consistent, we note that from the definition of the 
comoving curvature perturbation = —4' -|- Tiv and from equation (2.7) we have 

CL = X + 2x-b. (2.10) 

Thus, a long-wavelength solution is a gauge mode if on those large-seales the comoving curva¬ 
ture perturbation is constant in time. This means that the long-wavelength limit of a physical 
solution of Bardeen’s potentials is equivalent to a coordinate transformation and it is what 
Weinberg calls adiabatic modes [31]. Indeed, since ( is conserved outside the sound hori¬ 
zon which is nearly zero in a ACDM universe, this would hold up to scales where baryonic 
physics start being important. Note that this constancy of f must be satisfied only by the 
long-wavelength mode for our arguments to be valid; the short-wavelength modes can be at 
any scale. 


3 Observed galaxy overdensity in the presence of long mode 


In this section we compute the effect of a large-scale mode approximated by a constant plus 
a gradient as in Eq. (2.2). We start by defining the galaxy number density, we then include 
the effect of the long mode and we show that if the long mode has a characteristic scale 
that is much larger than all other scales in the problem (including the distance between the 
observer and the galaxies), its effect is zero as expected. 


3.1 Galaxy number density 


In galaxy surveys, the observed galaxy number density at a given redshift and angle on the 
sky is obtained by counting the number of galaxies, dNg within the observed volume. 


dV = dXldz 


nf\z,h) 


dAg(z,h) 

dH 


(3.1) 


From here, one defines the galaxy number overdensity Ag as 


nf\z, h) = ng(2;)(l -7 Ag(2;, h)). 


(3.2) 
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where hg(z) = (n°*’®(z)) is the observed number density of galaxies averaged over the angle 
at a fixed redshift. 

In order to study the effect of a long-wavelength mode on the galaxy number density, 
we must first relate it to the physical density. In the inhomogenous universe, the observed 
redshift and the propagation direction of the photons differ from the true redshift and angle. 
Therefore the observed volume defined in terms of the observed redshift and solid angle differs 
from the physical volume occupied by the source 

diphys = dV'obs(l + 5V) (3.3) 

Thus the observed number density of galaxies which is obtained by counting the number of 
galaxies within the observed volume differs from the physical number density by a factor due 
to the volume distortions 


nf\z,h) = nf^\z,h){l + 8V) . (3.4) 

Moreover, galaxy surveys are flux limited, i.e., only galaxies above some threshold luminosity 
are observed. This threshold luminosity is inferred from the observed flux 

L = 4Tr'Df{z)fohs, (3.5) 

where the luminosity distance T^i{z) is that in a homogenous universe and L is the inferred 
luminosity. However the physical luminosity of the galaxy is different than the inferred 
luminosity since the luminosity distance to the galaxy is not that of the homogenous universe 
and the photon propagation is affected by the fluctuations along the path 

^phys ^ ^T^ViizfU, = 47rPK2)Vobs(l + mf = L{1 + 5Vif (3.6) 

Therefore the observed number of galaxies at observed redshift z, angle n and above a 
threshold L is related to the physical number count as 

nf>%z, h, L) = nfy%z, h, L(1 + + 6V) 

= nf^^iz,h){l + 6V){l + t5Vi) (3.7) 

where t = 2 —■ In writing the second line in Eq. (3.7), as explained in the beginning 
of this section, we made two assumptions: first we only kept the linear contributions to the lu¬ 
minosity distance and volume perturbations. This is because if we neglect the strong lensing, 
the perturbations in the luminosity distance and volume are proportional to the gravita¬ 
tional potential and its gradient and thus their higher order contributions are subdominant 
with respect to the higher order contribution of the matter density. Thus we neglect these 
terms. Second in writing the contribution from the magnification bias, we assumed that the 
perturbations and the background have the same dependance on the luminosity. 
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3.2 Long-short couplings 

We now wish to compute the galaxy over-density measured in a small patch of the sky and 
in a small range of redshifts, in the presence of a mode of Bardeen’s potential which is much 
larger than the size of this patch. Following the arguments of section 2, this mode should 
have no effect on the local physics of that small region since it is equivalent to a coordinate 
transformation. However, we take here this large-scale mode to have a characteristic size 
which is smaller than the distance between the observer and the small patch, which means 
that it will indeed have an effect on the correlations measured by the observer. This effect 
is purely due to the fact that we measure the galaxy correlation functions in terms of the 
redshift and angular position of the galaxies rather than the physical separation between 
galaxies and their proper time. For this reason this non-zero correlation is a “projection 
effect”. 

In the presence of the long mode, the observed redshift and angle of the galaxies are 
modified with respect to those in the absence of long mode. First lets compute the effect of 
the large-scale mode on the observed angular position of the galaxy 

f,L = ^ =, (3.8) 

where = (5)^ is the projection of the photon 4-momentum orthogonal 

to the observer’s 4-velocity vP. At linear order the solution to the photon geodesic gives 

no = ne + Vq - f dr' V-^ 4'(r')) , (3.9) 

J TO 

where the direction of emission is defined as % = —PpTg)/^y^6jPP^^p)PP^p) and the 
superscript denotes a projection orthogonal to hg, e.g. u*-*- = {6^^ — h^h^)uT The second 
term on the right hand side corresponds to an aberration due to the motion of the observer 
and the last term is the lensing of the photon trajectory. 

The presence of the large-scale mode will contribute to the lensing integral in Eq. (3.9), 
and it will also cause the photon to be emitted in a different direction (in the coordinate 
system with the long mode). This change in the direction of emission can be computed 
simply by remembering that at the point of emission the large-scale mode is equivalent to a 
coordinate transformation which will transform the vector P^{t£) giving 

ns^ns = n£ - ^“^(r^) -F 2b-'-(T£ - to) , (3.10) 

where we used the fact that at zeroth order xg = — to)- Writing ^ and b in terms of 

the potentials as in Eqs. (2.8) and (2.10) and isolating the contribution of the long mode to 
the lensing integral in Eq. (3.9) we get 

no = no-VQ + {l+'H{T£-To))v-t{x£)-{T£-To)V^^Lix£)+ f dr'V-^ ($L(r') + ^l(t')) • 

J TO 

(3.11) 
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The redshift is defined as the ratio of the energies of a photon at emission and observation 

^ P^{rs)us,i 
P'^{to)uov ’ 

and a calculation analogous to the one performed for the direction gives 


z = z + {l + z) 


t-TE 

'HvL{xs)-'Hovo + ^L{xs)-^o-n-{vL{xs)-vo)- / dr'($'^ + T'^)(r') 

Jto 

(3.12) 

The observed number density of galaxies in the presence of the long mode is therefore given 
by 


nf\z, n, L) = nf^\z + Az, n + Ah)(l + 6VP^){1 + 

L 

— ^gi^) [1 + ^ (^1 + eAz + 6VP^ + , (3.13) 


9 d In 

where^ e = —gj—, and Ag has contributions from the luminosity distance and volume 
perturbations due to the short modes as well as the density perturbations due to the short 
modes. Note that Az = z — z and An = n — h are the corrections to the observed angle and 
redshift due to the long mode. Expanding Eq. (3.13) linearly the galaxy number overdensity 
is finally given by 


^g('2^) n)\L = Ag{z, h) + (1 + Ag(2;, h)) (^eAz + 

d d f \ 

(raj 

For our purposes it suffices to use the perturbative expressions for the quantities and 

at first order [5, 6, 32] 


JCW = - 2(4., + 5.,) + i'I'i + (^ + T) J., 

+ + (3.15) 


and 


= 


rU 


rro 


- 1 -VL ■h + ^L+ dr' (T'^ + 


(Tl + $l) (r) - 


1 

fro 

2 r — r(r') 

+ “ 

/ dr' 

-/VAo 

2 

Jr 

r rr[T) 


(3.16) 


^In the literature the evolution bias is often defined as /evo = 


d ln(a 


ndr£ 


therefore /evo = 3 — e/a. 
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where all the quantities without an argument are evaluated at the point of emission. 

The last step we need to take is to correlate the expression (3.14) with another short- 
scale Ag and a long-scale one, and go to harmonic space to obtain the consistency relation 
of the galaxy three-point function in the squeezed limit 


^Ag (zi,hi)Ag(z2,h2)Ag(2;3,h3)^ = (^Ag{zi,hi)d{z2,h2)'^ + (^A^{zi,hi)Az{z2,n2)'^-^ 


d 

+ (^Ag{zi,hi)An{z2,n2)'^ ■ (^Ag{z2,n2)Ag(z3,h3) 


+ ( 203 ), 


(3.17) 


where we have introduced 

d = {eAz + (3.18) 

Let us close this section by commenting on the corrections to Eq. (3.14). In general, we 
also expect corrections going like the curvature induced by the long-wavelength mode, that 
is {^LAg/y^^TH?'). Parametrically, these are subdominant with respect to the lensing and 
redshift space distortions which behave as (<l>iAg/AiA5?7^), but they are larger than the 
other relativistic corrections. However as we will discuss in section 4.2, there are additional 
suppressions for the lensing and the redshift space distortions. In spite of this, our result is 
a consistency relation in the sense that if the evolution of the universe is adiabatic (single- 
field inflation and a subsequent evolution that conserves ^ on the scales of interest) our 
relation captures the behavior of the squeezed limit as we take l/A^ to zero up to first 
order in Ag/A^. Indeed, the terms we computed could be mistaken for a non-zero local 
non-Gaussianity, which corresponds to a physical coupling between and <!>£:, or a signal of 
violation of the Equivalence Principle if they are not taken into account. Eurthermore, there 
are several approaches to compute these curvature terms in the literature, which correspond 
to computing the evolution of the short modes in a cnrved universe [33-35] (see also [36, 37]) 
but they can only be approximate since it is in general impossible to analytically compute 
the effect of curvature on the evolution of these non-perturbative short modes. 

Note that in real observations, one counts the number of galaxies within a finite redshift 
bin which have two effects. Eirst, the characteristic scale of the short mode is determined 
both by its angnlar scale and the assumed width of the redshift slices. Therefore the curvature 
terms are only negligible if \l ^ 5z/T-L. Second as we will explain in more details in section 
4, in order to relate the theoretically calcnlated bispectrum to the one measured in a snrvey, 
one needs to integrate over the window functions describing the redshift bins. Note that the 
consistency relation does not break down when performing this integration as long as the 
characteristic spatial scale of the redshift bin is smaller than the characteristic scale of the 
long-wavelength mode. This is easily achieved especially for large redshifts z^l. 

3.3 Validity checks 

Let us now perform various validity checks on our results. 
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• We expect a large mode to have no observable effect when its scale is larger than 
all the scales in the problem, i.e. the distance between the source galaxies and the 
distance between them and the observer. This is similar to what happens for the CMB 
[27, 29]. One way to check that our calculation is consistent with this fact is to perform 
a coordinate transformation that induces such a large-scale mode and check that the 
galaxy number over-density remains invariant. Let us begin by checking that this is 
indeed the case for ho, setting the observer’s position to zero xq = 0, we get 

no = no + i'^iro) - + 2 {ts-T o)b-^ + f dr' {t') - 2b) , (3.19) 

all the terms of the transformation cancel. An analogous argument holds for the 
redshift z, and we get 


z = 


z + {l + z) 


e'{T£) + X£ ■ tins) + h-^{ ts) 


e'{To) -h-i{ to) 



(3.20) 


where all the terms in the parenthesis cancel after writing the partial time derivative in 
terms of a total derivative in the integral dr = d/dr -|- h • V. For these cancellations to 
happen it is important to replace the Bardeen’s potentials and T inside the integral 
by an exact constant plus a gradient, Eq. (2.2), which does not hold if the long mode 
oscillates between the observer and the source, i.e. if the characteristic scale of the 
long mode is smaller than the distance between the observer and the source. We thus 
obtain 

rf^''{z,h)\L = h^{z) {I + l:^^{z,h)) + 2t5v\^^'^ . (3.21) 

Finally, the average in this equation hg assumes that the observer has access to the 
large-scale mode, such that the terms in the second parenthesis give zero when averaged 
over angles. However, when the mode is much larger than all the scales in the problem, 
the angular averages of these terms do not vanish. 


{nf^{z,h))\L = {h^{z)) (l + + {2t5v\^^)) , (3.22) 


but they are zero-modes on the observed patch and they can be reabsorbed in the 
average of the galaxy number density 


nf\z,h)\L = {nf\z,h))\L (1 + Ag(z,h)) . 


(3.23) 


We thus conclude that a mode with a characteristic scale much larger than all the 
scales in the problem has no effect on the Ag. We remark once more that when the 
large-scale mode is smaller than the distance between the source and the observer, these 
contributions will all be non-zero. 
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• As a further check of the validity of our calculation we take the super-horizon limit in 
equation (3.14) and work to first order in the long modes and zeroth order in the short 
modes. We expect to obtain the same results as the linear calculation of Ref. [5-8], 
which in a pressureless medium is given by 

= b6 m,sync -h (-efcsync + tSVi + SV) , 

where 5m,sync and 5^sync are the matter over density and redshift perturbations in 
the synchronous gauge. Indeed, Az = —6zsync and at super-horizon scales 5m,sync is 
suppressed by and we recover our result. 

• As remarked in the literature, the expressions for and are gauge invariant. 

The reader might worry that Az and Ah are not since they explicitly involve the 
velocity which is not gauge invariant. However, it is straightforward to write them at 
first order in terms of explicitly gauge invariant variables 


Ah = 


+ {rE- To) ) ^^Cl{x£) + ] 


'^('Tf) 


^{xe) 


r'^e 


+ dr'V^(4>L(r') + 4'L(T')) , (3.24) 


JTO 

Az = (l+z) Cl{xe)+^L{ xE)+^L{xE)-^fi-V{CL+^ l){xe)- [ dr' (4>'i(r')-FT'i(r')) 

^ J TQ 

(3.25) 

where we have dropped velocities and potentials evaluated at the observer as we will 
do in the rest of this paper. 


4 Spherical Harmonic Decomposition 

In both N-body simulations and data from galaxy surveys, the three-point function is often 
measured and analyzed in Fourier-space. This is a convenient decomposition if the galaxies 
have similar redshifts. However here we are interested in large-scales for the long mode where 
the relativistic effects can be important. Hence a single Fourier mode of large wavelength may 
include galaxies with significantly different redshifts. Due to the evolution of the universe, 
there is no translation invariance in the radial direction of the redshift space, i.e. z. Therefore 
the Fourier space bispectrum is not proportional to the Dirac delta: triangles do not close. 
We thus choose to decompose the three-point function in spherical harmonics, which could 
prove advantageous or even necessary for this type of analysis, 

(a£imi(5;i)a^2m2(^2)a£3m3(23))„ „ = dDidD2dD3y/^^yhi)y/2^2(h2)T/3m3(ra3) 

\ / J 

X (|Ag(zi,hi)Ag( 2 : 2 ,h 2 )Ag( 2 ; 3 ,h 3 )y (4.1) 

This basis naturally fits the sphere and does not have the problem of triangles not closing. We 
now compute each of the terms contributing to the above equation in turn. The contribution 
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of the terms in the first line of Eq. (3.17) are easy to compute and give 


('2^l)Q-£2m2 ('^2)®£3m3 (^ 3 ) ) « „ ^ ^ 




X +(2 0 3) , 


(4.2) 


where the Gaunt integral is defined as 

<3T,l7r = I Ye^mAn) Ye,mM Ye.mAn) 

^ (h h i 3 \ (il h h \ / (24 + I)(2^2 + 1)(24 + 1) 

\rni m 2 ms) \0 0 0 J \ dvr 


The third line of Eq. (3.17) is slightly more complex 


(ae,rm{zi)ai^m2{z2)aesm3{z3))^ E/^^^(hi)y/2m2(^2)T/3^3(h3) 

\ / Ei<^T2/3 J 

d d 

X ^(|Ag(zi,hi)I(z2,n2)) • ^^Ag(2;2,h2)Ag(z3,h3)^ + (2 o 3) (4.4) 

where 

I{z2,n2) ^^^1+ (-'H{t 2)] VL+ f -^^($L + TL)(r',h2r(T')). (4.5) 

\r(T2) ) Jro r{T') 

This can be further simplified by decomposing the two-point functions in spherical harmonics, 
giving 


^iimi (zi)ci£2m2 {z2)cLigm3 (-^-S, 


^l<SC^2/3 


D 


X Z2)C^t^'(^2, ^ 3 ) + (2 o 3). (4.6) 


The integral in the first line can be rewritten as 


1 

2 


dfl 


Aq (y;,^, (h)y;3^3 (h)) - (A^y/^^^ («)) y;3m3 («) - Ye^m, («) i^nYel^^ (n)) 
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(4.7) 


The second and third terms in the parenthesis give simply 

+1) + hih +1)). (4.8) 

We write the product of spherical harmonics of the first term as a Clebsch-Gordan decom¬ 
position 
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Putting everything together for the observed reduced bispectrum of galaxies, &£i,£ 2,^3 -^ 2 , ^ 3 ), 

defined as 

(^l)®t 2 m 2 ('2'2)Q.£3m3 (- 23 )^ ~ Q ^ (^ 1 ; ^ 2 ; -2^3)) (4-10) 

we obtain the main result of this paper, namely the consistency relation in the squeezed limit 
in multipole space 


, Ijm, &^i,£2,4(^1>^2,^3) = 




+ ^ (£i(4 + 1) - ^ 2(^2 + 1) + mz + 1)) Cff{z^,Z2) 


c;^"'^Hz2,^3) + (2f^3). 




J4JI) 


As pointed out before, since in real observation, one counts the number of galaxies 
within a finite redshift bin, in order to compare with observed galaxy bispectrum, one needs 
to integrate our theoretically calculated bispectrum in (4.11) over the three window functions 
describing the redshift bins 


^^,^ 2 , 4 (^ 1 ’^ 2 , 23 ) = j dz[ dz 2 d4 W{zi, z[)W{z 2 , Z 2 )W{z 3 , z'^) Z 2 , z'^) . (4.12) 

Given our expression in (4.11), it is straightforward to compute ^^{zi^Z 2 ,z^). 

4.1 The distant observer approximation 

In order to build some intuition, let us rewrite our expression (4.11) in the distant observer 
approximation where we choose to decompose in Fourier modes only the direction perpen¬ 
dicular to the line of sight. 


Ag(zi,/i)Ag(z2,/2)Ag(2;3,4)4 , „ „ = (2vr)^<5D(G +-^ + ^) 




C^^\zl,Z2,^l) 


d 


+ , Z 2 , G ) 7 ^ + 4 • (^1 , ^2 , 4 ) 

OZ2 


C^-^4z2,^3,4) + (2o3). (4.13) 


Note that the third term in the square brackets, when added to the corresponding one coming 
from (2 -H- 3), is 0(£f) in the equal time limit as expected due to momentum conservation, 
while the other terms are different from zero in that limit. Note also that the Newtonian 
consistency relations take into account only the redshift space distortions induced by the 
long-wavelength velocity, which are included here in Az. This Newtonian contribution is 
not zero in the equal-time limit. The same would hold if one also decomposes the direction 
parallel to the line of sight in Fourier modes. This is due to the fact that one performs 
measurements on the light-cone, and since there is no translation invariance in the radial 
direction in redshift space, the two momenta connected to this radial direction would not 
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Figure 1. Reduced observed bispectrum in the squeezed limit at equal-time (z=0.55) on linear scales. 
Note that the consistency relation is also valid when the short modes € 2)^3 o-fe in the non-linear regime. 
The vertical line indicates the linear scale up to which the numerical code is valid. 

cancel. However, this non-cancellation will be governed by how much the power spectrum 
changes with redshift; when it is a good approximation to ignore the evolution of the power 
spectrum this Newtonian consistency relation is expected to be zero as was previously found. 
Indeed, in the above expression it is proportional to a derivative of the power spectrum with 
respect to redshift. Let us mention, however, that the contribution from this term is very 
small due to the fact that velocity and density are not correlated in harmonic space. 

In the rest of this paper, we use the general expression given in Eq. (4.II) without 
resorting to the flat sky approximation. 

4.2 Second-order limit 

As an illustration, we plot in figure 1 different terms contributing to the reduced bispectrum 
appearing in the relation (4.11) in the case where the short modes are linear. The second- 
order computation of the three-point correlator of the observed galaxies should reproduce 
this result in the squeezed limit. We do this using a modified version of the code CLASSgal 
[38]. The result is presented in the simple case of no galaxy bias b = 1, no magnihcation bias 
t = 0 and no galaxy evolution /evo = 0 for equal redshifts of 0.55. The vertical line denotes 
the linear scale at z = 0.55 of £ ~ 83 up to which the numerical code is valid. 

One could be surprised by the small size of the cross correlation between Ag and Az. 
In the Newtonian approximation and working at a hxed time slice, the terms proportional 
to Az are indeed believed to be the dominant effect. However, the term is suppressed in 
harmonic space due to the fact that the Bessel functions corresponding to the Ag term and 
the h ■ vp term are ji and j'^ which oscillate at a similar frequency but out of phase and 
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Figure 2. Effective induced by the relativistic effects computed in this paper. The three different 
lines correspond to different magnification bias parameters. We have assumed e = 3/a for all cases. 
This estimate is valid for z > 1. At lower redshifts, integration over the redshift bins can change this 
result. 


therefore cancel. In Fourier space, this can be intuitively understood from the fact that 
the regions where the velocities are at the maximum correspond to the minima (in absolute 
value) of the density and vice versa. Moreover the lensing contribution is smaller than one 
would expect parametrically. This is due to the fact that the integration along the line of 
sight tends to average out the perturbations. 

5 Effective local non-Gaussianity parameter from non-linear general rela¬ 
tivity corrections 

As mentioned in previous sections, the terms proportional to ‘hAg would induce a signal in the 
observations that might be misinterpreted as primordial non-Gaussianity with a local shape. 
In order to see this, let us estimate the squeezed limit of the galaxy bispectrum induced by 
the primordial local non-Gaussianity. From a Newtonian calculation and assuming linear 
biasing we have^ 

' ' h-l ?h3 

where /nl ~ 1 is the non-linear non-Gaussianity parameter parametrizing the strength of the 
primordial non-Gaussianity with a local shape. To see the impact of the non-linear relativistic 
contributions, we compute the ratio of Eq. (4.11) to where the potentially highly 

similar expression is obtained in Ref. [23] for the dark matter over-density bispectrum and it is expected 
to hold at the non-perturbative level. However, it is not clear whether this is also true for galaxies. 
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non-linear short-mode power spectrum cancels such that our result is fully non-perturbative'^. 
We then isolate the terms for which this ratio is independent of i.e. those terms in Eq. (4.11) 
with the same t'-dependence as Eq. (5.1). In order to compare with the observed bispectrum 
in a galaxy survey, one needs to integrate over the window functions defining the redshift bins 
before taking this ratio. Unless the window functions are delta functions, which for a realistic 
survey they are not, this integration will change both the amplitude and the Udependance of 
the bispectrum. However, due to the nearly constant behaviour of the transfer functions at 
z > 1, the non-linear power spectrum of the short modes still cancels in the ratio. Therefore 
this integration does not modify our estimate of for those redshifts. Eor lower redshifts, 
one can explicitly carry out the integration, though the results would be only perturbative in 
the short modes. It should be noted that since the signal for the primordial non-Gaussianity 
is strongest at higher redshifts, removing the GR effects that mimic this signal is of crucial 
importance at those redshifts. 

The exact value of will depend on redshift and the magnification bias. We plot its 
value as a function of redshift in figure 2 for three different values of the magnification bias. 
Erom these results we deduce that at redshift z = l.b and for a vanishing magnification bias, 
the value of the effective is about —3.0. A qualitatively similar result was obtained in [5] 
for the power spectrum. They found a very similar dependance on redshift and magnification 
bias as well as a value for of the same order of magnitude. Since they considered a 
different observable, the exact values are different though their behavior is qualitatively the 
same. 

Also, these relativistic effects could also induce terms that would be mistaken for a 
violation of the Equivalence Principle in a naive analysis [39, 40]. 

6 Conclusions 

In this paper we have exploited the fact that a long-wavelength mode of Bardeen’s potential 
is an adiabatic mode to write a non-perturbative relativistic relation for the squeezed limit 
of the galaxy number overdensity three-point function. Adiabatic modes are physical large- 
scale solutions to the dynamics that can be cancelled by a residual gauge transformation. 
The latter is realized non-linearly by Bardeen’s potential. Thus, computing their effect on 
short-scale physics reduces trivially to a change of frame. 

The effect of these transformations on the observed galaxy number density cor¬ 
responds to a change in the redshift and the direction of the photons at observation which 
in turn changes the relation between the observed volume and the physical volume, and a 
magnification bias. When computing the galaxy number overdensity Ag, one should also 

^This cancelation does not happen for the term proportional to /dz. However, as one can 

easily check at second-order, this term has a different ^-dependence than Eq. (5.1). This is because the 
derivative with respect to redshift acting on the short mode power spectrum, acts on the Bessel function and 
will induce a different ^-dependence. A similar ^dependence will also arise at higher orders 
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include the change of the average hg as a result of distortions to the redshift, the so-called 
evolution bias. 

Curvature and tidal contributions, proportional to second derivatives of the gravita¬ 
tional potential, are not captured by our argument and are expected to be larger inside the 
horizon than some of the relativistic corrections which we keep, even though they would 
naively be expected to be parametrically subdominant with respect to the redshift space 
distortions and lensing induced by the long mode. However, redshift space distortions, which 
were also obtained in a Newtonian computation, are small due to the requirement that the 
small-scale modes be on similar redshifts (though they will be non-zero even at equal red¬ 
shift), and the lensing is suppressed since it averages in the trajectory of the photon from 
emission to observation. We argue that the relativistic corrections we compute can however 
be distinguished by their dependence on the large scale as one takes the squeezed limit, 
and they serve as a consistency check for our cosmological model. Indeed, the violation of 
the Equivalence Principle at large scales or the presence of primordial local non-Gaussianity 
would induce a squeezed limit with the shape of these contributions. They need to be kept 
in order not to confuse them with these deviations from the vanilla cosmological scenario. 
The exact value of the effective local non-Gaussianity parameter from these non-linear GR 
corrections depends on the redshift and the magnification bias. At redshift of z = 1.5 and in 
the absence of magnihcation bias we estimate = —3.0. 

Finally, since the universe evolves and there is no translation invariance in the radial 
direction of redshift space, Fourier-space triangles do not close and we thus present our 
expressions after decomposing them in spherical harmonics. Notably, we find that the New¬ 
tonian redshift space distortions induce a term that does not go to zero at equal redshifts, 
but is however very small. 

There are many directions in which this work can be extended. It would be interesting 
to see to which level of accuracy planned large scale structure surveys could put constraints 
on these consistency relations. The prospect is promising since one can use arbitrarily small 
scales for the short modes given that the consistency relations are non-perturbative. However, 
the effect of systematic observational effects must be taken into account. Another interesting 
direction would be to perform a similar calculation for the galaxy two-point function including 
contributions from the dark matter bispectrum, thus affecting the measured bias and its scale 
dependence which is sensitive to the squeezed limit; such computation would include all the 
relativistic and non-linear effects. 
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